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Abstract: The fiber length trait (FLT) of 538 individuals from nature birch population in Maorshan region, Heilongjang, China were meas¬ 
ured, of which 100 individuals were selected as representative variety of correlated fragments screening with random amplified polymor¬ 
phism DNA (RAPD) technique. In total of 20 RAPD primers were tested through multiple regression analysis between amplified strip and 
the character behaviors, and a correlative segment BFLR-16 was obtained. The correlation coefficient between BFLI-16 and FLT was 
0.6144, with the significant level of 1%. BFLI-16 was then cloned, sequenced and transformed into SCAR marker. The percentage of iden¬ 
tifying long fiber birches by this SCAR was more than 92. The result indicates that the SCAR markers has high specificity for the long fiber 
individuals and is highly linked with the gene controlling the character of fiber length, and its existence is significantly correlative with the 
increase in the fiber length. 
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Introduction 

Betula platyphylla (silver birch) has the characteristics of well 
updated, rapid growth, strong adaptability, wide distributing, 
high content of cellulose with low content of lignin and white 
color etc. (Li et al. 1995). It is recognized as a kind of su¬ 
perordinary and short-period material for paper pulp. Studies 
showed that the fiber length was an important index for the me¬ 
chanical strength of paper, which affects paper’s qualities. There¬ 
fore, one of important aims for researching the fiber length trait 
(FTL) of birch is to select new species of birch. Molecular label¬ 
ing technique has an incomparable advantages compared with 
traditional breeding methods, and it has been widely used to 
establish genetic linkage maps of trees and quantitative trait lo¬ 
cus (QTLs) (Cervera et al. 2001; Wu et al. 2000). However, the 
QTLs of a special trait selected from a single-population is 
greatly limited due to the restriction of labeling theories and 
technical methods. And the comparability between the results of 
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the QTLs of different labs selected from different populations is 
low (Zhang et al. 2003; Grattapaglia et al. 1994; Jermstad et al. 
2003). The researched populations and analysis methods of 
QTLs need improvement. The main random molecular labeling 
techniques used in the research of QTLs such as RAPD (Random 
amplification polymorphism DNA), ISSR (Inter-simple sequence 
repeat) and AFLP (Amplified fragment length polymorphism) 
have the disadvantages of low specificity and stability, and the 
domain amplified is unknown. Therefore, the gained QTLs can 
hardly be used in trees’ labeling aided to selecting breeding. In 
the SCAR labeling technique, a pair of specific primers is de¬ 
signed according to a known sequence of the genome. Then a 
restriction region of the genome is amplified with high specific¬ 
ity, which is used as a specific labeling to the genome. This 
method overcomes the limitation of the RAPD markers tech¬ 
nique and is widely used for locating relative genes of a correla¬ 
tive traits and establishing genetic maps (Meletto 1996; Garcia 
1996; Ohmori 1996). 

In this study, a fragment closely correlated with fiber length 
trail of birch was selected using the random amplified polymor¬ 
phism DNA (RAPD) technique, then sequenced and cloned, 
finally transformed into a stable SCAR marker. The aim of this 
study is 1) to provide a tool for high accuracy marker assistant 
selection in longer fiber birch, and 2) to set a good background 
and a new method for marker assistant breeding on other species. 

Materials and methods 

Plant materials 

A total of 583 individuals of silver birch were randomly selected 
from natural birch forest in Maorshan Forest Experiment Farm of 
Northeast Forestry University, Heilongjiang, China. The core 
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was collected with an awl at diameter at breast height (DBH) 
from south to north. Tender leaves were collected with scissors at 
the same time. The leaves and the core were ensured to be from 
the same tree in order to make markers. The collected leaves 
were kept immediately in the refrigerator at -20°C. To identify 
the validity of the segments, the verified materials of 200 birch 
individuals selected randomly were collected at the Harbin Ex¬ 
perimental Forestry Farm of Northeast Forestry University with 
the same methods. These 200 individuals did not have any ge¬ 
netic relationship with the natural birch forest in Maorshan. 

Methods 

The determination of fiber length 

The pith was cut into small pieces with the same length and 
thickness as matchsticks. Then it was kept in 30%-50% nitric 
acid with slight potassium chloride solution. This can decompose 
the pith into cellulose. The fiber length was measured under 
optical microscope with an enlarged multiple of 40 to 70. In each 
sample, 10 pieces of fiber selected randomly were measured and 
the average value standed for the firer length of this sample. 

Genomic DNA extracting 

The improved CTAB methods were used to extract the whole 
DNA of the silver birch (Zhan et al. 2005). 

Analysis of RAPD 

The RAPD reaction system was determined through orthogonal 
design methods, on which the best reaction system of silver birch 
was gained. The 20-pL reaction system includes 2.4-pL Mg 2+ (25 
mmol/L), 0.8-pL dNTPs (2.5 mmol/L), 1.5-pL primer (50 
pmol/pL), 0.1 -pL Taq DNA polymerase (5 U/pL), and 2.5-pL 
template DNA (15 ng/pL). Amplification reaction was performed 
on a thermal cycler (MJ-9600PCR), programmed as follows: 3 
min at 94°C for initial denaturizing, then 38 cycles of 94°C for 
30 s, 38°C for 30 s, 72°C for 90 s, and extended at 72°C for 7 
min, finally keeping it at 4°C. 

The individuals with the longest and the shortest fiber length 
were taken as samples and 20 RAPD primers were screened for 
good polymorphism. 100 genomic DNA of the silver birch was 
used as template for the amplification reaction, and PCR prod¬ 
ucts were separated by electrophoresis at 150 V for 2 h, through 
the agar gel of 1.5%. After taking pictures in the UVP gel imag¬ 
ing system, the diverse strip was noted, 1 stands for having strip 
and 0 for none. 

The calculation of correlation coefficient 

Taking the genotypic value of signed segments as the independ¬ 
ent variable and the fiber length of 100 selected birch individuals 
as variable, we establish the multiple regression equation as fol¬ 
lows: 

y n = V + b n x n + b i2 x i2 +k + ... b m x m 

Where, y n is the fiber length of birch individual (n = 1,2.100), 

ju the regression constant, b ni the correlation coefficient of 


signed segment with fiber length (/ is the sequence number of the 

signed segments, i = 1,2,3,4,5,6,7.), and X ni represents the 

genotypic value of a signed segment of the silver birch having 
the sequence number of n (X, 0/1). 

The relativity of fiber length with molecular markers is a re¬ 
spective multiple regression coefficient. When the regression 
coefficient is positive, it indicates the existence of signed seg¬ 
ments is positively correlated with the increase in fiber length. 
When negative it indicates the existence of signed segments is 
negatively correlated with the decrease in fiber length. The sig¬ 
nificance level of the regression coefficient represents the degree 
of gene controlling the character of fiber length linked with the 
signed segments. If the significance level is high, it indicates the 
degree of gene linked is tight. Using the software package SPSS 
12.0, we analyzed the regression coefficient of all signed seg¬ 
ments, from which we selected the signed segments with high 
significance level. Those segments were reclaimed and trans¬ 
formed to SCAR markers. 

Reclaiming and sequencing the segments with a high signifi¬ 
cance level 

PCR products were separated by electrophoresis on the agar gel 
of 1.5%. The strip notably correlative with fiber length was cut 
from the gel. Reclaimed it with the reclaiming reagent boxes 
(purchased from bao biological company of dalian). It was 
cloned to the carrier of PMD18-T and transformed to the TOP 10 
strain. After checking with PCR, it was sent to be sequenced by 
bao biological company of dalian. 

Transformation and verification of SCAR markers 

Suitable primers were designed with software package primer 
5.0 according to the results of sequencing. The primers were 
synthesized by shanghai biological company. The 20-pL ampli¬ 
fied system includes 1.2-pL Mg 2+ (25 mmol/L), 0.4-pL dNTPs 
(2.5 mmol/L), 0.4-pL cis-primer (50 pmol/pL), 0.4-pL 
trans-primer (50 pmol/pL), 1-pL Taq DNA polymerase (1 U/pL), 
1-pL template DNA (40 ng/pL), and 2-pL 10 x PCR buffer. Am¬ 
plification reaction was performed on a thermal cycler 
(MJ-9600PCR) programmed as follows: 4 min at 94°C for initial 
denaturizing, then 35 cycles of 94°C for 45 s, 56°C for 45 s,72°C 
for 80 s, and extended at 72°C for 7 min, finally keeping it at 
4°C. PCR products were separated by electrophoresis at 120 V 
for about 2 h, through the agar gel of 2%. And the pictures were 
taken in the UVP gel imaging system. 

Results and analysis 

Selection of birch individuals 

The fiber length of 583 individuals was analyzed with the soft¬ 
ware package of SPSS. Results showed that the average fiber 
length of the 583 individual was 4.73, which exhibited a typical 
normal distribution (Fig. 1). The results were in consistence with 
the feature of phenotype distribution of quantity genes control- 
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ling characters. If the fiber length is larger than the mean value 
(4.73), then a plus standard deviation (0.36) is given; oppositely, 
minus standard deviation is given. Base on the rule, we respec¬ 
tively chose 50 individuals with the longest and shortest fiber 
length to do the genome’s discrepancy analysis. The fiber length 
of 200 individuals selected from the Harbin Experimental Forest 
Farm of Northeast Forestry University was analyzed with the 



Compaiatrze length 

Fig. 1 Distribution of liber length of 583 birches 


software package of SPSS. Results showed that the average fiber 
length was 5.46 and the population showed a normal distribution 
(Fig. 2). The results were consistent with the feature of pheno¬ 
type distribution of quantity gene controlling characters. Based 
on the above rule (the mean value (5.46), the standard deviation 
(0.60)), we selected 25 individuals with the longest fiber length 
and 25 individuals with the shortest as the verified population. 



f ig. 2 Distribution of liber length of 100 birches 


Screening of relative RAPD signed segments 

A total of 20 RAPD primers were used to analyze the genome of 
100 birch individuals selected. Each primer had 0 to 1 diverse 


strip. The primer R-16 (TGCCGAGCTG) gave one different 
strip and the correlation coefficient was 0.6144. This segment 
was named “BFLR-16” with a molecular weight of 1200 bp or so 
(Fig. 3). 


Ml 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 



Fig. 3 Electrophoresis of amplilied products willi primer R-16 in some long fiber individuals M DL2000 
DNA molecular marker. 1-25: Long liber indiv idnals; the arrow indicates the correlated fragment 


The reclamation and sequence of the correlation segments 

The correlative segment BFLR-16 was reclaimed. The amplified 
reaction was conducted again with this segment as template and 
R-16 as primer. The molecular weight of the amplified segment 
was about 1200 bp (Fig. 4), matched with the objective segment, 
therefore, can be used to conjunction and clone directly. The 
reclaimed segments were connected to the carrier of PMD18-T 
and then transformed to the TOP 10 strain. After cultured in the 
plate, 3 to 4 single white colonies were selected randomly to be 
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cultured in the rocking bed. The plasmids were extracted. The 
amplified reaction was conducted using R-16 as primer to detect 
the correctness of the inserted segments (Fig. 5). Results indi¬ 
cated the recombinant plasmid had only one strip after electro¬ 
phoresis and the molecular weight was consistent with correla¬ 
tive RAPD segments with long fiber trait. It demonstrates the 
RAPD specific segment has been correctively inserted into the 
plasmid. The bacterium including RAPD specific segment cor¬ 
relative with long fiber trait was sent to bao biological company 
to be sequenced. The sequence was as follow: 












Journal of Forestry Research (2008) 19(4):288—292 


291 


TGCCGAGCTGCTCAAATCCATTAATCCCAATCCATGATAATGGGTAAATTGGGTAACAATTGCATTGTTTAAAAAAACTCCC 


AAAATTAACACGAAGTGTCGTCGTGGTGCAAACTCCCAAATGATAATTAAGGTGAACTGTCATTTGTACCAACAAATGTT 


TCTAGGGAAAAAGAAAATAAACACATAAATTCACGACTTAATAAGTAATATCAGAATAAATACAACAGAATCAATCTAAGTT 
A AGTA ATATTT CTTT CT C A A A A AT C A ACTT C ATTAG A A A A A A A ACTATTTACCTT G ATA ACTTTAC A AGTTA A AGCCTCC AG 
A AT GC ATTA AC AT CTA A A ATATACTATT CTAGC ATA A ATAC A AC AGTACCC A A A A A ATT CAT GTA A AT CT GCT CTT C AGG ACT 
TTAGCCAGCAACTCAAAACCTCACATTAACTCTAATTTGCAATTTGCTAATATGTAAAACAATAAATTGCCAACGAAATTTA 
GCCCATGCAGTGAGCTCGGTTCGC ATAAATACAAC AGTACCC ATGTTGG AAAC AAA ATATCATTAAGGCTTAGTTTACAGTT 
CAT GTA A ACTA ATTT GCTA ATAT GTA ATT CTTATTAT G AG A A A A ACTAG AGT GGG A AT G A A ATTTAT C AGTT GTA A AGCTTATT 
TCTTCACCACACCTAAATTTTATCCAGAATGCATTAACACCCAAAAAATATCTAAAATATACTATTCTAAAAACCCCGAGGCT 
AAATTCATGTAAATCTGCTCTTCAGGACTTTAGCCAGCAACTCAAAACCTCACATTAACTCTAATTTGCAACCCCGAGGCTA 
ATTTGCTAATATGTAAAACAATAATGTTTCTAGAAAATAAACACATAAATTCACGACTTAAAAGCCTCCAGAATGCATTAACA 


T CTA A A ATATACTATT CTAGC ATA A ATAC A A AC A ATA AT GTTT CTAGGG A A A A AG A A A ATA A AT CTGCTCTT G AC C AC A ATA 


GGCAGCTTCACCGAGGCTAGAGTGGGAATGAAATTAGTAACCTCACATTAACTCTAATATCAGTTGTAAAGCTTATTCACA 


CCTA A ATTAC ATTA ACT CTAT CTT C ACC AC ACCTA A ATTACT CTA ATTAC ATTATTAT GGTATT C A AGTA ACCT C AC ATTA ACT C 
TAATTACATTAACTCTATGCAACCTTATCAGCTCGGCA 


M 


M 




250 

100 


I ig. 4 Electrophoresis of amplification products with 
primer R-16 for callback M, DL2000 DNA molecular 
marker 1-4. Electrophoresis of amplification products 
for callback 


Fig. 5 Electrophoresis of amplification products with primer 
R-16 for recombines TOP 10 M, DL2000 DNA molecular 
marker 1-4, Electrophoresis of amplification products for 
recombines TOP 10 
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Fig. 6 Electrophoresis of the amplified products with primer BFL i n 50 individuals birch Note; M. 
1)1 2000 DNA molecular marker: 3-25, long fiber individuals, 26-50, short fiber individuals: The 
arrow indicates the correlated fragment 


The verification of SCAR marker 


According to the sequenced results, the specific primers with 20 


bp were selected with the software of primer premier 5.0. This 
RAPD marker was transformed to SCAR marker. The sequence 
of BFLR-16 transformed to SCAR-R16 was as follow: BFLR-3F; 
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5 'GGTGAACTGTC ATTTGTACC35 TGAAGCTGCCTATT- 
GTGGTC'. The length of the amplified segment was 440 bp. The 
genome of 50 birch individuals selected was used as template to 
be conducted the PCR with the designed primer. PCR products 
were separated by electrophoresis on the agar gel of 2%. There 
were 23 out of 25 individuals with longer fiber length having 
diverse strip (Fig. 6). The probability of long fiber individuals 
amplified was 92%. One out of 25 short fiber individuals was 
amplified. Therefore, it indicated the SCAR markers had high 
specificity for the long fiber individuals and may be highly 
linked with the gene controlling the character of fiber length. 

Discussion 

The forest has many important economic traits such as tree 
height, diameter at breast height, leaf area, and fiber length that 
are all quantitative characters (Zhang et al. 2002). These charac¬ 
ters are controlled by many chromosomal loci (QTLs), usually 
distributed in different provenance sources or family constella¬ 
tions. The segregation population used to locate QTLs comes 
from the cross of two inbred lines, but this kind of segregation 
population is hardly to be obtained for perennial trees. Even if 
we can use the FI progenies to substitute the segregation popula¬ 
tion, only part of the controlled characters can be identified 
(Zeng et al. 2006). This research gained 583 birch individuals 
with unknown genetic background. After doing multiple correla¬ 
tion analysis about the genotype and phenotypic characters of 
100 typical individuals, we obtained the segments correlative 
with long fiber. Not only did we avoid the restriction of obtain¬ 
ing the QTLs hard for a single-population, but also there was 
hope to reveal all the QTLs systematically correlative with 
birch’s long fiber trait as many character genes were gathered in 
the natural birch forest with complex genetic background. 

The analytical method of correlation coefficient can be used to 
analyze and identify the sites of quantitative characters (Xu et al. 
1994). The molecular marker correlative with the coefficient 
obtained by the method was consistent the results of variance 
analysis (Xu and Zou 2002). Compared with the mapping soft¬ 
ware MAPMAKER/QTL based on the maximum likelihood 
statistical single interval graphic method, the analytical method 
of correlation coefficient can be used to do correlation analysis 
about any solitary molecular marker. Using this method, there is 
no need to do several molecular markers to establish detailed and 
complete genetic maps and the material of birch can be both a 
specific segregation population and a combination of randomly 
selected individuals. This is helpful for analyzing systematically 
the expression state of the quantitative characters. Therefore, the 
method connecting the character of birch’s fiber length with 
corresponding genomic marker is feasible. The markers corre¬ 
lated with birch’s long fiber trait were obtained by multiple re¬ 
gression analysis. Also, it is possible to breakthrough the situa¬ 
tion of retarded progression on marker aided selective breeding. 

The RAPD method is easy to operate, but it is random marker 
and has the lowest rate of repeatability. Therefore, it’s necessary 
to transform the obtained RAPD signed segments to the stable 
SCAR markers. In this research, we successfully transformed a 
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RAPD segment “BFLR-16” to SCAR marker with an identifica¬ 
tion efficiency of 92% for long fiber birch, which can be used to 
aid selective breeding of long fiber birch. But the identification 
efficiency of SCAR marker is not absolute (Tang et al. 2005; Shi 
et al. 2001). Therefore, there is still need to make systematically 
research about every step of the process RAPD segments trans¬ 
forming to SCAR. 
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